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Hadamard spectroscopy has earlier been used to speed-up multi-dimensional NMR ex- 
periments. In this work we speed-up the two-dimensional quantum computing scheme, by 
using Hadamard spectroscopy in the indirect dimension, resulting in a scheme which is faster 
and requires the Fourier transformation only in the direct dimension. Two and three qubit 
quantum gates are implemented with an extra observer qubit. We also use one-dimensional 
Hadamard spectroscopy for binary information storage by spatial encoding and implemen- 
tation of a parallel search algorithm. 



I. INTRODUCTION 

The use of quantum systems for information processing was first introduced by Benioff In 
1985 Deutsch described quantum computers which exploit the superposition of multi particle states, 
thereby achieving massive parallelism P]. Researchers have also studied the possibility of solving 
certain types of problems more efficiently than can be done on conventional computers ^, 0, Q] ■ 
These theoretical possibilities have generated significant interest for experimental realization of 
quantum computers . Several techniques are being exploited for quantum computing and 
quantum information processing including nuclear magnetic resonance (NMR) ^, . 

NMR has played a leading role for the practical demonstration of quantum gates and algorithms 
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. In NMR, individual spins having different Larmor frequencies and weakly coupled to 



each other are treated as individual qubits. The unitary operators needed for the implementation, 
have mostly been realized using spin selective as well as transition selective radio frequency pulses 
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and coupling evolution, utilizing spin-spin (J) or dipolar couplings among the spins |l3l.ll4l.ll5l.ll€|. 
The final step of any quantum computation is the read out of the output. In NMR the read out is 
obtained by selectively detecting the magnetization of each spin or by tomography of full density 
matrix [l^. [l^. It was first proposed by Ernst and co-workers, that a two dimensional experiment 
can be used to correlate, input and output states, which is advantageous from spectroscopic view 
point In two-dimensional quantum computation (2D QC), an extra qubit (observer qubit) 

is used, whose spectral lines indicate the quantum states of the work qubits The 2D spec- 

trum of the observer qubit, gives input-output correlation of the computation performed on the 
work qubits 19]. The 2D spectrum is therefore more informative than a one dimensional (ID) 
spectrum. For example, in ID NMR QIP, the spectrum after the SWAP operation, performed on 
the equilibrium state of a homonuclear system, is identical to the equilibrium spectrum. However 
the same operation, performed using a 2D experiment, contains the signature of SWAP gate [2^. 
The observer qubit can also be used to prepare a pair of pseudo pure states 2ll |. 



logic gates and several algorithms are implemented by 2D NMR 
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The quantum 



20|,l211. Recently, 2D NMR 



las also been used to address the decoherence free sub spaces, for quantum information processing 



22|. 



Multi-dimensional NMR spectroscopy is often time consuming, since each indirect dimension 
has to be incremented to span the whole frequency range, and the desired digital resolution 2^ . 
Several experimental protocols have been developed to accelerate the recording of multi-dimensional 
spectra. These include, single scan experiments in the presence of large gradients, GFT, Covariance 

"^nnn 



spectroscopy and Hadamard spectroscopy 
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27l . l28l . I29I . |2d] • The Hadamard spectroscopy, 



proposed by Kupce and Freeman, has the advantage that one can simultaneously label, various 
transitions of the spectrum by applying a multi-frequency pulse j^^, ^ suitable decoding 

followed by a Fourier transform only in the direct dimension yields a 2D spectrum [27 Jj^ L This 
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leads to a large saving in time for experiments having a small number of transitions 
this paper we use Hadamard spectroscopy to speed-up the two dimensional quantum computing 
scheme ll 



Information storage and retrieval at the atomic and molecular level has been an active area 



of research |33, 
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39 1 . Khitrin et.al demonstrated that, the spectrum of 



dipolar coupled spin cluster can be used to store large amounts of information, which can be used 

]. Alternately, it 



for photography and implementation of parallel search algorithm 
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has been demonstrated that, spatial encoding under a linear field gradient can also be used for 



above purposes 
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4J]. In this work, the one-dimensional Hadamard spectroscopy 



3l| | is used 



under spatial encoding, to store the information and to implement a parallel search algorithm. The 
proposed method has the advantage that, once the Hadamard encoded data is recorded, one can 
write any binary information array (sentence), and search for any code or alphabet in that array. 
The main emphasis of this paper is to demonstrate the use of Hadamard encoding in the field of 
NMR information processing. 

In section (II), we outline the Hadamard method for 2D-NMR QIP along with the conventional 
method. In section (III) we implement, various 2D-gates on 3 and 4-qubit systems. In section (IV), 
we implement parallel search algorithm by using Hadamard spectroscopy under spatial encoding 
and section (V) contains the conclusions. 

II. THEORY 



Quantum computing using two-dimensional NMR can be described by transformations in the 
Liouville space. For a spin- 1/2 nucleus, having two orthogonal states |0) and |1), the longitudinal 



polarization operators can be written as 



IS 



, 1 0\ /O 0\ 1 1/10 

Jo = |0)(0| = , /i = |1)(1| = ], and h = -(/o - h) 



a) 



A product state \ip)= |001...0) of a coupled spin-1/2 nuclei, can be represented in the Liouville 
space by a density matrix a, obtained by the direct product of longitudinal operators. 



Inlnh ....In. 



'0 — -'O-'O-'l 



(2) 



In 2D-NMR QIP 19], an extra qubit (observer qubit) is used, whose transitions represent the 
quantum states of the work qubits (computation qubits). Thus an (N+l)-qubit system can be used 
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for N-qubit computation, treating the zeroth qubit as the observer qubit. The thermal equihbrium 
state of observer spin I*^, can be represented in the Liouville space as [l^ . 

< = + {Illl-I^) + + (3) 

where the superscript indicates the qubit number, with the observer qubit represented by the letter 
'O'. 

In the following we describe the conventional and the Hadamard 2D methods (Fig^, for a three 
qubit system, under the NOT operation on both the work qubits, during the computation period. 
The schematic energy level diagram of a three qubit system and the spectrum of the observer qubit 
are given in Fig EI The transitions of the observer qubit, which represent the quantum states of 
the other two qubits (work qubits), are labeled as |00), |01), |10) and A NOT operation 

performed during the computation period, interchanges the states |0) and |1) of both the work 
qubits. 



(A) The Conventional Method 

As shown in Fig.(^^), the observer spin is first allowed to evolve for a time ti during which the 
work qubits remain in their initial states, after the frequency labeling period ti, the computation 
is performed on the work qubits, followed by the detection in t2 period. A two-dimensional Fourier 
transform gives the 2D spectrum of the observer qubit, which shows the input-output correlation 
of the computation, performed on the work qubits. 

For a three qubit system, the equilibrium state of the observer qubit I*^, can be written as, 

= i?[{ihil) + {ihiD + {iliD + {lllDl (4) 
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The pulse sequence given in fig^^, transforms cr^ as, 



(7r/2)^ 
> 


X 


[{iiiD + mil) + + {IIP,)] 


ti 

> 


X 


[cOs{uJQQti){llll) + COs{uQiti){llll) + COs{L0iQt{){llll) + COs{L0iiti){llll)] 


> 




[cos{uJooti){I^li) + cos{ujoiti){l^lf) + cos(wioti)(/|/o) + cos{ujnti){lllf)] 


Unot 
> 




[cos{u;ooti){lllf) + cos{iuoih){lll^) + cos{LOioh){I^lf) + cos{LOuh){I^li)] (5) 


{n/2)y 
> 


X 


[cos{uJooh){lllf) + cos{ujoiti){lll^) + cos{uJioti){I^lf) + cos{uJuti){Ioli)] 


;> 


X 


[cos{uJooti)cos{u;nt2){lllf) + cos{uJoiti)cos{uJiot2){lllo) 



+ COs{LJioti)cOs{LJoit2)iIolf) + COs{LJiiti)cOs{uJoot2)iIolQ)] 



Fourier transform performed in both dimensions on the above signal, gives a two dimensional 
spectrum, where input and output states are given in Fi and F2 dimensions respectively. The 
time consuming part of this method is the large number of ti increments, needed to achieve the 
required spectral width and sufficient resolution in the Fi dimension. Quadrature detection in the 
Fi dimension, further doubles the number of experiments. 

(B) The Hadamard Method 

In this method (Fig03), the sequence (7r/2)-ti-(7r/2) - Gz of FigdA., is replaced by a multi- 
frequency vr (MF-vr) pulse on the observer qubit. Instead of ti increments of method (A), the 
pulse sequence of Fig^B, is repeated k times, where k = 2^ , is the number of transitions of the 
observer qubit. In each of the k experiments, the multi-frequency vr pulse is differently encoded, 
according to the rows of a k-dimensional Hadamard matrix. For a two work qubit case (k=4), four 
experiments are performed with four different encodings of the vr pulse, given by the four rows of 
the four dimensional Hadamard matrix (Fig. IHK), where '-' and '-|-' in the matrix corresponds to 

'tt pulse' and 'no vr pulse' respectively. For example, H h - means, the vr pulse is applied only on 

2nd and 4th transitions of the observer qubit. The output of the four experiments (FigEt'), under 
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the NOT operation on both the work qubits, can be calculated as follows, 
Experiment (1): 

< I? mil) + m) + m) + m)] 

S^^I^ lO[cos{iUut){lllf) + cos{iUwt){lll^) + COs{LUoit){I^lf) + cos{iuoot){l^l^)]- 
Experiment (2): 

^^^^ I? mii) - m) + m) - m)] 

^ lO[{llll) - (Illl) + {lllD - illll)] (7) 

lO[cos{ui^t)[l\ll) - cos{uiot)[l\ll) + cos{u^it)[llll) - cos{u^ot){llll)]; 
Experiment (3): 



miiiii) + {iiiD - m) - [iiiD] 

lO[cos{u;nt){lll!) + cosMilll^) - cos{uJoit){I^lf) - cos{uJoot){I^ll)]; 
Experiment (4): 



lO[cos{u:iit){llll) - cos{uJwt){llll) - cos{iOoit){I^lf) + cos{iOoot){Ioli)]; 



where (tt)'*-'^'!^'"^ means, a vr pulse is applied on and \lm) transitions of the observer qubit. 
Each of the four experiments (eqn.s 6 to 9) generates a composite response of the computation, 
performed on the work qubits. However, the different encoding pattern applied in each experiment, 
provides a decoding method of extracting the output state, individually for each of the input 
states. The decoding is obtained by the transpose of the Hadamard matrix. The decoding of the 
experiments 1, 2, 3 and 4, for the input states (Fig|33)> lOO), [01), [10) and |11) are respectively 
given by. 
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(1) + (2) + (3) + (4), 

(1) - (2) + (3) - (4), (10) 
(1) + (2) - (3) - (4), 
(1) - (2) - (3) + (4). 

A two-dimensional spectrum of the computation can be constructed, by inserting the decoded 
data (time domaiii) at suitable frequencies in the Fi dimension following a Fourier transform in 
the F2 dimension [23|. For example, in the present case, the decoded data (eqn.lO) for the input 
states |00), |01), |10) and |11) are inserted in the Fi dimension, respectively at the frequencies ujqq, 
<^oi, ^10 and wii (Fig. 121), following a Fourier transform in the F2 dimension, yielding the desired 
2D spectrum. 

It is to be noted that, the Hadamard encoding can also be achieved by J-evolution. For example 
in the above case (eqns. 6 to 9), the observer qubit can be represented in terms of product operators 
(respectively for eqn.s 6 to 9, figEK) as. If, (If /^), (/f li), and (if li/f), after the Hadamard 
encodin g (M F-tt pulse) . Each of these product operators can also be prepared by using J-evolution 



method 



While the conventional method (Fig^A.) needs a minimum number of ti increments for a sat- 
isfactory resolution in the Fi dimension, the Hadamard method (Fig^3) inherently yields high 
resolution in Fi dimension and needs only a small number of experiments, equal to number of 
transitions of the observer qubit. It may be noted that for an N work qubits system, the number 
of transitions of the observer qubit (for weakly coupled spins with all resolved transitions) is 2^, 
thus for small number of qubits (up to about 9 qubits) the Hadamard method is advantageous. 
The Hadamard method can also be used for 2D implementation of quantum algorithms 2(1. l21 1. 
It may be added that the Hadamard method does not change the scaling of quantum computing 
nor does it change the scaling of any algorithm. 
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III. EXPERIMENTAL IMPLEMENTATION OF QUANTUM GATES 

(a) Two qubit gates 

The system chosen for implementatfon of two qubit gates, is C2-F3/, where the three fluorines 
can be treated as three qubits. The fluorine spectra are shown in Fig. ^ The transitions of the 
observer qubit are labeled as |00), |01), |10) and |11). 

The two-qubit N0T(1,2) gate (NOT on qubits 1 and 2) is implemented (Fig. Et) using method 
(A), with 128 ti increments, a recycle delay of 20 seconds (~5Ti) and 2 scans for each increment, 
resulting in a total experimental time of 126 minutes. The two-qubit N0T(1,2) and several other 
gates are implemented by method (B), shown in Figl^j. The Hadamard encoding, shown in FigEJi, 
is achieved by MF-vr pulses of duration 100 ms.. The unitary operators and computation pulses 
for various two-qubit gates are given in l2dl]. The NOP gate is a unit matrix, hence the output 
states are same as input states, N0T(1,2) interchanges |0) and |1) of both the work qubits, SWAP 
gate interchanges the states |01) and jlO), and CNOT(l) gate interchanges the states |01) and |11). 
Each 2D gate of FigJSJa, is recorded in four experiments, which for same recycle delay as in method 
(A), takes the total experimental time of less than 2 minutes. 

(b) Three qubit gates 

The four fluorines of 2-amino, 3,4,5,6-tetra fluoro benzoic acid, can be used as four qubits. The 
one-dimensional spectra of observer qubit {I^) and work qubits (/^, and I^) are given in (FigEI, 

Since for this system the {I^) spin has 8 transitions the Hadamard method (Fig^B) for 2D 
QC gates requires 8 experiments as outlined in FiglSb- Due to small separation of the frequencies, 
FigEl ( 5 Hz, as compared to 40 Hz in Fig. 0]), the MF-vr pulse needs about 600 ms.. Hence the 
Hadamard encoding, in this case, is achieved by using J-evolution method explained below. 

The magnetization of the observer qubit, after the Hadamard encoding (Fig^a), can be repre- 
sented as, , /f /f, /], I^njl /f, /f/l/f, /i/f, (FigEb). The pulse sequence 

(V2)^-(l/2Jo,)-(vr/2)0 is used to prepare I^Il, where the evolution is with respect to Joi- The 



9 



product operator I^Illi {i / j) is prepared by the pulse sequence (7r/2)^-(l/2Joi)-(l/2Joj)- 
(7r/2)^, and I^lljlll is prepared by the pulse sequence, (7r/2)^-(l/2Joi)-(l/2Jo2)-(l/2Jo3)- 
(7r/2)5^. Each 1/2 Joi evolution is achieved by applying selective vr pulses simultaneously on the 
observer and the i^^ qubit, in the middle of the evolution period (l/2Joi). 

A NOP gate is implemented, which requires no operation during the computation, is shown in 
FiglZK- NOT(l) and N0T(2) are implemented by applying a vr pulse, respectively on and 
(FigEbic). The Toffoli gate, which is a universal gate for reversible computation, is achieved by 
applying two transition selective vr pulses on transitions |0110)-|0111) and |1110)-|1111) (I7|l). Each 
2D gate shown in Fig[71 takes about 2 minutes, for a recycle delay of 8 seconds. The conventional 
method with 128 ti increments, takes about 60 minutes (spectrum not shown). 



IV. PARALLEL SEARCH ALGORITHM 

Information storage by NMR, was suggested almost 50 years ago by Anderson et.al [3j]. This in- 
volves the excitation of various slices of the isotropic liquid (for exarnple H2O), under the z-gradient, 
using series of weak radio- frequency pulses followed by spin echo [Sj]. Khitrin et.al demonstrated 
that, the multi-frequency excitation of dipolar coupled spectrum of liquid crystal, enables a par- 
allel (simultaneous) storage of the information, at the atomic level 

Q- In ref.QQ. 

it is shown 

that, one can imprint the information written in a binary code, where '0' and '1', in the frequency 
space corresponds to no excitation and excitation respectively. A 215 bit sentence is written, where 

each alphabet is assigned a five bit string, for example a=l (00001), z=26(11100) and blank 

space as 00000 3. It is further shown that, one can perform a parallel search on the binary 
information array (sentence) using six bit-shifted multi-frequency pulses, to search for a letter, in 
a string of letters [43] . 

Recently, it has been demonstrated that, spatial encoding can also be used for information 
storage and parallel search using the single resonance of a liquid such as H2O in the presence of a 



linear field gradient 



42 



4J]. Spatial encoding involves radio- frequency (rf) excitation at multiple 



frequencies in the presence of a linear magnetic field gradient along the z-direction. Spatial encoding 
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gqrithm using XOR operation 



4J]. The XOR search requires 



was also used by Sersa et.al, to excite an arbitrary three dimensional patterns, using x, y, and z 

gradients Q, Q, Q- 

NMR photography and a parallel search al; 
have been implemented by spatial encoding under z-gradient 
only two experiments, in which the first experiment is used to record the sentence and the second to 
record the ancilla pattern of the letter to be searched. The absolute intensity difference spectrum 
of these two experiments is obtained, followed by the integration of the peak intensities of the 
5 bits corresponding to each of the letters. The pattern of the integrated intensities yields the 
zero intensity, only at those places where the required letter is present. The XOR search not 
only searches the required letter, but also searches the letter having the complementary code. For 
example, as shown in [4^, the letters "o = 01111" and it's complementary "p=10000" can be 
searched simultaneously. 

In this work, the spatial encoding is conjugated with Hadamard spectroscopy [s^l- We synthesize 
256 phase encoded multi- frequency 7r/2 pulses, each of which consists of 256 harmonics and excite 
the 256 slices of the water sample, under z-gradient (figlHl)- The phase encoding of the harmonics, 
in each of the pulses, is given by the rows of the 256 dimensional Hadamard matrix, where '+' and 
'-' in the matrix, corresponds to the phases y and -y respectively. Thus the application of the pulse 
under the gradient (figEl)) creates either or -Ix magnetization of each slice. The 256 pulses are 
used to record the Hadamard encoded data of 256 slices, by using the pulse sequence of figlHI As 
seen in the previous sections, the Hadamard encoded data can be suitably decoded to obtain any 
element (frequency) of the Hadamard matrix. It should be noted that, in this case the Hadamard 
encoding is performed on the Ix magnetization of various slices, whereas in the 2D gates (Fig^D, 
section (II)), the encoding is done on the z magnetization of the observer qubit transitions. The 
Hadamard encoded data of 256 slices, stored in 256 separate files, can be suitably decoded to write 
any binary information, which requires a maximum of 256 bits. 

The Hadamard encoded data is suitably decoded, to write the sentence, "the quick brown fox 
jumps over the lazy dog" (Fig|3\,B). The XOR search is performed to search a letter "u", the 
ancilla pattern for the letter "u" is decoded in Fig.ljUp), and the difference spectrum of Figsl^V 
andEP, is shown in Figlj^I)). Integration of absolute intensities of the peaks of each of the letters 
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in Fig.(jUp), are shown in Fig.(j^). The zero intensity in Fig. ((35), indicates the presence of letter 
"u". Maximum intensity is observed for the letter "j", which has a code complementary to letter 
"u". The advantage of method is seen in the next example, in which the same data is used to 
write another sentence and search code for another letter. Only the decoding pattern is different 
in example of Figs llUt V. contains another sentence, "principles of nuclear magnetic resonance", 
which consists of 200 bits of information. The letter "e" is searched by using the XOR search. 



The encoded data can also be used for 256 bits of NMR photography 



encoding 




3J, |40|, |4lj, 143] . The spatial 



4J] has the advantage that the relaxation of all the slices (bits) is uniform. It may be 



pointed out that in the liquid crystal method [43] as well as J-coupled systems |46| all the lines are 
not independent and perturbation of one line can cause disturbance in other lines of the spectrum, 
which forms the basis of the Z-COSY experiment j^]. On the other hand, in spatial encoding 
method the pattern is inhomogeneous broadened and parts of the spectrum can be independently 



perturbed |4! 



V. CONCLUSIONS 



In this paper we demonstrate the use of Hadamard encoding for (i) two-dimensional quantum 
information processing and (ii) for parallel search using spatial encoding. For (i), this method con- 
verts the 2D experiment to a small number of ID experiments requiring the Fourier transformation, 
only in the direct dimension. The required encoding can be achieved by using multi-frequency tt 
pulses or J-evolution method. For (ii) the Multi-frequency excitation and detection of the water 
sample, in the presence of z-gradient, maps the magnetization of various slices, to the frequency 
space. Each slice is treated as a classical bit which can exists either in or 1, which in the fre- 
quency space respectively correspond to no excitation and excitation. The Hadamard encoded data 
is suitably decoded for the information storage and implementation of parallel search algorithms. 
It will be interesting to use the Hadamard spectroscopy, for information storage using x,y and z 
gradients. 
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FIGURE CAPTIONS 

(1) .(a) Pulse sequence for 2D NMR QIP. is the observer qubit, and /^,..../^ are work 
qubits. During ti period the input states of the work qubits are labeled followed by the compu- 
tation, and signal acquisition during the t2 period. A two dimensional Fourier transform results 
the 2D spectrum of the observer qubit, where the input and output states are given in Fi and F2 
dimensions respectively, (b) pulse sequence for 2D Hadamard NMR QIP, k experiments are per- 
formed, where k is the number of transitions of the observer qubit. In each of the k-1 experiments, 
the MF-vr pulse is applied on k/2 transitions (explained in text and figEI)- The results of the k 
experiments can be suitably decoded, to obtain the output state of the computation, individually 
for each of the input states. A two-dimensional spectrum is obtained by inserting the decoded data 
at suitable frequencies in the Fi dimension, following the Fourier transform in the F2 dimension. 

(2) (a) Schematic energy level diagram of a three qubit system, and deviation populations of 
the equilibrium state, (b) The equilibrium spectrum of the observer qubit, whose transitions are 
labeled as the quantum states of other two qubits. 

(3) (a) and (b) are Hadamard matrices, which are used to implement two and three qubit 
gates respectively (fig. EtjIZI), by using 2D Hadamard QIP (figQ^). Each of the columns of the 
Hadamard matrix are assigned to the transitions of the observer qubit. In the matrix '-|-' and '-' 
corresponds to no pulse and vr pulse respectively. The product operators, associated with each of 
the encodings, are also given in the last column. 

(4) Fluorine spectra of C2-F3/. The three fluorines form a three qubit system. The chemical 
shifts of work qubits with respect to observer qubit, are i7i=11807 Hz, ^2= -17114 Hz, and the 
J-couplings are Joi= 68.1 Hz, Jo2= -128.8 Hz, and Ji2= 48.9 Hz. The transitions of the observer 
qubit l'^ represent the quantum states of the work qubits (/^ and /^). The relative signs were 
determined by selective spin tickling experiments [2^. 

(5) (a) Implementation of two qubit N0T(1,2) gate, on a three qubit system (figEJ,by using 
conventional method given in fig^A.; 128 ti increments are used, with 2 scans for each increment and 
a recycle delay of 20 sec, resulting in a total experimental time of 126 minutes, (b) Implementation 
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of two qubit quantum gates by using 2D Hadamard QIP (fig^). Each 2D gate is recorded in four 
experiments, taking the total experimental time of less than 2 minutes. The encoding of the M.F. 
TT pulses, in each of the four experiments, is given in fig|2K- NOP gate requires no pulse during the 
computation. N0T(1,2) is implemented by applying selective vr pulses on both the work qubits 
and I^. Swap gate requires six transition selective vr pulses on transitions |110)-|111), |010)-|011), 
|101)-|111), |001)-|011), |110)-|111) and |010)-|011), CNOT(l) requires two selective vr pulses on 
transitions |001)-|011) and |101)-|111). The phases of the vr pulses during the computation, are set 
as (x, -X, -y, -y), in order to reduce the distortions due to pulse imperfections. 

(6) . Fluorine spectrum of tetra-fluro benzene. The four fluorines form a four qubit system, 
where is the observer qubit, whose transitions are labeled as the quantum states of the three 
work qubits (I^, and I^). The chemical shifts of work qubits with respect to observer qubit, are 
Oi=13564.2 Hz, 172=6845.8 Hz, ^3 = -5261.2 Hz, and the J-couplings are Joi=10.5 HZ, Jo2=20.5 
HZ, Jo3=6 HZ, Ji2=9.5 HZ, Ji3=22.7 HZ and J23=21.9 HZ. 

(7) Implementation of three qubit gates by using 2D Hadamard QIP (figO^)! the ID spectra of 
the observer qubit /*~^ and the work qubits I^, and are shown in figEl Each 2D spectrum is 
recorded in 8 experiments. The Hadamard encoding, in each of the 8 experiments, is achieved by 
J-evolution method (explained in text). NOP gate is a unit matrix, hence the output states are 
same as the input states. NOT(l) and N0T(2), interchanges the states |0) and |1), of the I'st and 
2'nd qubits respectively. Toffoli gate interchanges the states |0) and |1) of the third qubit (I^), 
provided the other two work qubits (/^ and I^) are in state |1). 

(8) . Pulse sequence for the implementation of parallel search algorithm. The multi- frequency 
7r/2 pulse is obtained by modulating the Gaussian pulse with 256 harmonics and the phase modula- 
tion for each of the harmonics is y or -y. 256 multi- frequency 7r/2 pulses are synthesized, which differ 
from each other, only in the phase modulation, which is according to the rows of 256-dimension 
Hadamard matrix, where -|- and - in the matrix corresponds to the phases y and -y respectively. 
The duration of the M.F. pulse is 30 m.s., and the gradient strength is 25 Gauss/cm. The 256 ID 
spectra obtained individually from 256 pulses, are independently stored for suitable decoding, for 
desired information and search as shown in Fig. l9l and [Till 
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(9) . (A) spectrum obtained from the Hadamard decoding of the 256 experiments (hglHI), which 
represents the sentence (B), the quick brown fox jumps over the lazy dog. The sentence (B), 
consists of 215 bits or 43 ciphers, where each cipher is a five bit string, with a = 00001, b = 
00010,.. ..,z = 11100, and space = 00000. The "0" and "1" corresponds to "no excitation" and 
"excitation" respectively. (C) Ancilla pattern for letter u = 10101, "uuu....u (215 times). (D) The 
difference spectrum of (A) and (C). (E) The heights of the bars represent the integration of peak 
intensities of each of the five bit string of spectrum (D), the zero intensity (represented by arrow), 
indicates the presence of letter "u", and the intensity of letter "j" having a complimentary code 
01010 (represented by *), is 4.92 units (theoretically 5 units). This method is known as XOR 
search. 

(10) . The Hadamard encoded data is also used to write another sentence, "principles of nuclear 
magnetic resonance" (B). An XOR search is performed for letter e = 00101 (C), and the results are 
given in (D) and (E). The letter z having the complimentary code 11010 is absent in (B). Hence 
there is no line of intensity 5 units in (E). The next intensity is of 4 units, which correspond to 
letters "j", "r" and "x", whose code differ by 4 units from that of letter "e". However, only letter 
"r" is present in "B", occurring thrice and marked by * in (E). 
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